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.4.^ Abstract. This article is concerned with the representation of curves by 

Cw means of integral invariants. In contrast to the classical differential invari- 

H ants they have the advantage of being less sensitive with respect to noise. The 

I I integral invariant most common in use is the circular integral invariant. A ma- 
jor drawback of this curve descriptor, however, is the absence of any uniqueness 

fSJ result for this representation. This article serves as a contribution towards clos- 

^ ing this gap by showing that the circular integral invariant is locally unique in 

[^^ a neighbourhood of the circle. The proof is an application of Riesz-Schauder 

l/~^ theory and the implicit function theorem in a Banach space setting. 

. 1. Introduction 

f^ In many applications one faces the challenge to model objects, or parts of objects, 

^~^ in a mathematical framework. As an example, one important task is to extract 

\ I an object from a given data set and manipulate it in a post-processing step in 

^ order to obtain further information. Typical applications include medical imaging, 

L> object tracking in a sequence of images, but also object recognition, where the post- 

%^ processing step consists of the comparison of the extracted object with a database 

C^ of reference objects. Similarly, such a comparison can be necessary in medical 

imaging in order to distinguish between healthy and diseased organs. To that end, 

however, one has to be able to decide whether two given objects are similar or not. 

This requires a representation of the objects that makes the application of standard 

similarity measures possible. 

Finding a suitable representation of the object of interest, depending on the type 
of application, is crucial as a first step. For simplicity, it is often assumed that the 
object is a simply connected bounded domain, allowing for the identification of the 
domain with its boundary. From a mathematical point of view, this assumption 
reduces the complexity of the representation. In addition, there exists a larger 
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number of descriptions of boundaries than of domains, and, consequently, more 
mathematical tools to analyze the geometry of the underlying objects. 

In 2D a common approach is to encode the contour of an object by the curvature 
function of its boundary curve. This approach has, for instance, been used in 
[8], where the authors set up a shape space of planar curves, where the shapes 
are implicitly encoded by the curvature function. The main advantage of using a 
differential invariant — the most prominent representative being the curvature — 
to represent an object is the well investigated mathematical framework of this type 
of invariants (see [1, 10, 12]). 

Since all kinds of differential invariants are based on derivatives, they suffer from 
the shortcoming of being sensitive with respect to small perturbations. To bypass 
this shortcoming Manay et al. [11] proposed to use integral invariants instead of 
their differential counterparts (see also [6, 7, 16]). Integral invariants have similar 
invariance properties as differential invariants, but have proven to be considerably 
more robust with respect to noise. Their theory, however, is not that well investi- 
gated as opposed to the theory of their differential counterparts. 

Beside the classical approach of differential invariants and the novel approach of 
integral invariants, there exist several other concepts for encoding an object. For 
instance, in [3] the authors use the zero level set of a harmonic function, which is 
uniquely determined by prescribing two functions on the boundary of an annulus, 
to encode the boundary of a 2D object (see [4] for a generalization to compact 
surfaces in 3D). A similar encoding of the object by a function is given in the article 
of Sharon and Mumford [15]. Here, the authors first map the 2D object, which is 
supposed to be a smooth and simply closed curve, to the interior of the unit disc 
in the complex plane via the Riemann mapping theorem. This conformal mapping 
is composed with a second one, generated out of the exterior of the original object, 
and the composition is restricted to the boundary of the unit disc. Thus, the final 
mapping, which the authors call the fingerprint of the object, is a diffeomorphism 
from the unit circle onto itself. 

One of the challenges in object encoding is the question of uniqueness of the en- 
coding. More precisely, in many applications, e.g. object matching, the correspon- 
dence between the object and its encoding should be one-to-one. Thus, a thorough 
investigation of the operator that maps an object to its encoding is needed. In case 
of the encoding by a harmonic or conformal mapping — if possible — uniqueness is 
well known. Also for the encoding of an arc length parameterized curve by its cur- 
vature function, it is known that one obtains a one-to-one correspondence between 
the curve and its encoding (up to rigid body motions). One even has a complete 
characterization of the set of functions that arise as curvature functions of a class of 
sufficiently regular curves (see [2]). For integral invariants the situation is different; 
the cone area invariant, first introduced in [6], is an injective mapping independent 
of the space dimension, but its application is limited to star-shaped objects. In 
contrast, for the circular integral invariant, which is the integral invariant most 
common in use, there exists no proof for the uniqueness conjecture so far. 

This article is a contribution towards this goal: We first prove the Frechet dif- 
ferentiability of the integral invariant in a neighborhood of the circle, seen as a 
mapping from C*^'*"^ to C*^, fc > 0. Then we derive a local uniqueness result by 
proving the local injectivity of the Frechet differential in case fc > 1. The proofs of 
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these results are based on the impHcit function theorem on Banach spaces and an 
appHcation of Riesz-Schauder theory. 

2. Setting 

Let Emb be the space of all continuous embeddings from S^ to M^. Then every 
curve 7 G Emb has a unique interior, denoted by Int(7). Following [6, 11], this 
allows us to introduce the circular integral invariant: 

Definition 2.1. For given r > we define the circular integral invariant 

of a curve 7 G Emb as 

IrhK^p) ■■= area (5^(7 (</?)) n Int(7)) , 
where Br{p) denotes the ball of radius r centered at p G M^. 

The circular integral invariant behaves well under several group actions: 

• Ir is invariant with respect to Euclidean motions: For A G SE[2] we have 

Ir[Ao-f] =Irh] ■ 

• Ir is equivariant with respect to reparametrizations: For every homeomor- 
phism ^■. S^ —i' S^ we have 

I J. [7 o cf)] = Ir [7] o $ . 

• For every scalar i > we have 

Ir[t-f]^t^Ir/tb] . 

The observations above suggest to consider the integral invariant on the space 
C of all curves modulo Euclidean motions and reparametrizations. Moreover, we 
assume as an additional smoothness property that the considered curves are of class 
C'^, k > 1. Then it makes sense to use the following representation of C, as it avoids 
working with equivalence classes of curves. 

Definition 2.2. Denote by Ck C Emb, A; > 1, the space of all curves 7 G C'^(S'^; M^) 
satisfying the following conditions: 

• 7 has constant speed, i.e., there exists a constant Cj > such that ||7(((5)|| = 
Cj for all if & S^. 

• 7(0) = (1,0) and 7(0) = (0, c^), where we identify the circle S^ with the 
interval [0, 27r). 

• 7 is an embedding, that is, jiif) ^ lii^) for all Lp ^ ^. 

For the proof of our main theorem we need the following result from differential 
geometry concerning the manifold structure of Ck '■ 

Theorem 2.3. For fc > 1 the space Ck is a smooth submanifold of the Banach space 
of all C'^ -curves from S^ to M^. Its tangent space T^Ck at a curve 7 G Cfc consists 
of all C'^ -curves a with 

(ct((^), 7((^)) = C for some C eR, a{0) = (0, 0) and ((t(0), 7(0)) = . 
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Proof. The proof of the submanifold result is similar to [14, Thm. 2.2]. In our case 
the situation is less complicated, as we only deal with C'^-curves instead of Sobolev 
curves. The constant speed parameterization yields the condition 

2C = dMii^) + e&i'p),7i'p) + ea{^)) = 2{a{ip),^{ip)) . 

The remaining constraints follow directly from the initial conditions. D 

For k > 2 we obtain the following characterization of the tangent space T^Ck ■ 

Lemma 2.4. For k >2 the tangent space of Ck at a curve 7 € C^ with curvature 
function 

K^{^) .- -^ 

7 

consists of all C^ -curves a{tf) = a{Lp)'^{(f)^ + b{(f)'^{ip) satisfying: 

• 6((/3) = 6(0) + a{ip)n^{'.p)c^. 

• a(0) = 6(0) = 0. 

• d(0) = 0. 

Proof. Theorem 2.3 and the fact that (7((^),7((y3)) = imply that there exists a 
constant C S M such that 

C = (cr((^), 7((^)) = (7(^3), h{ip)^{ip)^ + a((p)7((^)-L + 6(<yj)7(v) + Kv)l{v)) 

= a(v)(7(<^),7('^)^) + Kv)c^^ = -a('^)4'«7('^) + Kv^^ ■ 

Using the initial conditions for a, we obtain the initial conditions for a and h and 
the value of C = 0^6(0). D 

We are now able to formulate the main result of this article: 



Theorem 2.5. The circular integral invariant Ir'. Ck+2 -^ C^{S^\M.), k > 1, is 
Frechet differentiahle on a neighborhood lA C Cfc+2 of the circle of radius R > r/2 
and its tangential mapping is injective on this neighborhood. 

In particular, it follows that Ir is locally unique on lA in the sense that each ^ £lA 
has a neighborhood V^ with Ir [7] 7^ Ir [7] for every 7 € V^ . 

Remark 1. The condition on r to be smaller than 2R is necessary, because oth- 
erwise the circular integral invariant in each point ip is constant equal to R^n, the 
area of the circle, and the same holds for any sufficiently small deformation of the 
circle which preserves the area. 

3. Frechet Differentiability of the Circular Integral Invariant 

In the following, we show that Ir is Frechet differentiahle and derive an analytic 
formula for Ir and its derivative /^ valid in a neighborhood of the circle. As a 
first step, we recall the following result on the differentiability of the composition 
mapping. 

Lemma 3.1. For every k > the composition mapping 

Comp: C''+2 X C'' ^ C'' , 

{f,9)^ f°9, 

is Frechet differentiable. 
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Proof. A Banach space valued mapping / : X — > y is called weakly differentiable, if 
the differential df : X x X ^ Y exists and is continuous. In contrast, it is Frechet 
dijferentiable, ii df : X — > L{X,Y) is in addition continuous for the operator norm 
on the image space. 

According to [5, 13], the mapping Comp: (7*^+^ x C'"' — !• C'' is weak C^. A 
straightforward calculation shows that a weak C'^-mapping is Frechet differentiable, 
which concludes the proof. D 

Remark 2. We will need two different types of derivatives for the formulation 
of our results: First, Frechet derivatives in the function space C'^, and, second, 
derivatives of functions / € C''{S^) with respect to their argument ip E S^. In 
order to highlight this difference, we use the following notation: For a function 
F : C^ ^^ C^ , we denote by F' : C'^ ^^ L{C^ , C) its Frechet derivative. In contrast, 
if / G C*^ , then / denotes its derivative in the parameter space. 

In order to make the notation less cumbersome, we omit the argument ip in 7(1/'), 
(y{(p) and similar expressions if the argument is clear from the context. 

Lemma 3.2. For each fc > there exists a neighborhood V C C'^+^(S'^; M^) of the 
constant speed parameterized circle of radius R > r/2 such that the following hold: 

(1) For each 7 G V and ip E S^ the circle Br{j{^)) intersects the curve 7 in 
exactly two points, denoted by ^[pj{ip)) and^{m^[ip)'). Herem^{(f) denotes 
the previous intersection parameter and p-y{ip) the next one (see Figure 1). 

(2) The mappings 

m: C''+^{S^;M.^)^V>m^{S^), p: C'=+2(5^M^) ^ Diff'^(S'i) , 

are Frechet differentiable. The derivatives in direction a G C'^+^(S'"'^; M^) 
are given by 

, ^ (g-g(p^),7(p-y)-7) ^, ^ (g-(7(m^),7(m^)-7) 

'' {i{p-i)ii{p-,)-i) ' '' {iim^),'y{mj)--f) 

Here Diff (S^) denotes the group of C -diffeomorphisms on the unit circle. 

Proof. Denote by 70 : -S*^ — > M^ the constant speed parameterized circle, that is, 
7o((/3) = (i?cos((^),i?sin((p)). Then, for every p G 5*-^, the circle Br{'yo{'f)) inter- 
sects 7o precisely at the two points (^Rcos{p ± ??), Rsm{(p ± ^)), where 

^ = arccos(l-^). 

Now define the mapping 

F: C''+^{S^;R'^) x Diff''(S'i) ^ C'=(S'^M) , 

(7,d)->|7(d)-7p-i?'- 

Obviously, the mappings pj and m^ we are searching for satisfy F{'y,p.y) = and 
F{'y,mj) = 0. In particular, the equation F{jQ,d) = has the two solutions 
Pjg{ip) :^ p + •& and mjg{ip) := ip — "d. Now, Lemma 3.1 implies that the mapping 
Comp: C"=+2(S'i; K2) ^ Diff'=(5'i) -^ C''{S^;R^), and consequently also F, is Frechet 
differentiable. Moreover, it is easy to see that the derivative of F at (7o,P7o) in 
direction (0, r) G C''+^{S^;R^) x C''{S^;R) is given as 

F'ilo,P'ro)i^'T) = 2(7o(p^J -7o,7o(P7o))^ = 2sin(?9)T, 
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which is obviously an isomorphism of C''{S^;M.), as the assumption R > r/2 > 
imphes that sin('!?) y^ 0. Thus the imphcit function theorem on Banach spaces 
(see [9, Sec. 1.5]) imphes the existence of a neighborhood V C C'^~^^ {S^ ■,M.'^) of 
7o and unique Frechet differentiable mappings ra, p:V ^> Diff (5*^) satisfying the 
equations i^(7,p-y) = = F{"f,m.y). The formula for the directional derivative of p 
at 7 in direction a now follows from the fact that 

= d^F{-f,p^){cr) = 2(7(p^) - 7,ct(p^) - a) + 2p'^{a){-f{p^) - 7,7(^7)) ' 

which is a simple application of the chain rule. The formula for m' (a) can be 
derived analogously. D 

Theorem 3.3. For each fc > there exists a neighborhood V C C'^"'"^(S'^;M^) of 
the circle of radius R > r/2 such that the following hold: 

(1) For each 7 G V the circular integral invariant Iri'j] can be written as 
(1) m (^) ^ I £ hW - 7, iW-)d^ + y arccos ( ^^^^^^ " ^;j(-^) ' ^^ 

(2) The circular integral invariant 

Ir- V^C''{S^-,M.), 

is Frechet differentiable. Its derivative in direction a E C''^'^{S^;M.'^) is 
given by 

2/; [7] (a) 

= 2 r'(a(V'),7W^)d^ 

- (c^. 7(^7)^ - li^i)^) + (7(^7) - 7, 0-^7)^) - (7(^7) - 7, cr(m^)^) 

-K7b7)-7,7(P7)") ^7^^-^;"^^-^-;^ 

(7(^7)' 7(^7) -7) 

frr- 

(7(771^)- 7, 7(m^) )- 



,±\ ( g-g(m^),7("^7)-7) 
(7(m^),7(TO^) -7) 



\2 



\/r^ - (7(^7) - 7, 7(^7) ~ 7) 

-(7b7)'7("i7)-7) 



(cr-cr(p^),7(p^)-7)/^^ 



(7(P7)'7(P7)-7) 

+ {<^{P-,) - cr, 7(^%) - 7) + (7(^7) - 7, cr(m^) - cr) 



(7b7)-7,7("i7)) 



_^ (g-o-(m^),7(w7)-7) 
(7(m^),7(m^)-7) 

Proof. Under the given assumptions for R, r and V, Formula (1) can be easily de- 
duced from Figure 1. A term by term investigation of Formula (1), using Lemma 3.1 
and Lemma 3.2, shows the Frechet differentiability of I^ on V. 
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7(to^(V3)) 




l{Pyi^)) 



Figure 1. Sketch of the derivation of the analytical formula for 
the circular integral invariant assuming two points of intersection. 

To calculate the differential of /^ we treat the two terms of Formula (1) separately. 
For the first term we obtain 



d^ 



(fT(^) - CT, 7(V;)^) + (7(V) - 7, aiyj)^)d^ 
+ iliPi) - 7, i{Pi)^)p'^y{<^) - (7(^7) - 7, 7(wT,)^)r7i^(cr) 



M^)-a,7W^>dV'- 



(7W,a(7/;)^)dV 



+ (7(^7) - 7, criPi)^) - ilimj) - 7, ^^(^7)^) 

+ (7(^7) - 7,7b7)^)P7(^) ~ (7(^7) -7,7("i7)^)TO^(cr) 



"^ 2(a(^), 7(^)^)rfV' - (^, 7(^7)^ - 7(m^)^) 

+ (7(^7) - 7, criP-r)^) - (7("^7) ~ 7, ^^(^7)^) 
+ iliPj) ~ J,J{Pi)^)p'-yi<^) ~ (7(^7) -7,7('7i7)^)TO^(cr) 
A simple application of the chain rule yields for the second term 

- a,^ arccos (^ (7(^7) - 7,7(^7) - 7) 
2 V r^ 



{'i{P-/)Pyi<^) +<^{P-,) - cr,jim,^) -7> + (7(^7) -7,7("^7)"^7('^) + ^^(^"7) ^ ^) 



^J'^-fK 



2r 2^r4 - (7(p^) -7,7(771^) -7)2 
Using the formulas for the intersection parameters, we obtain the desired result. D 

In the special case where 7 equals the unit circle the lemma above reduces to: 

Lemma 3.4. Let 7 e C'^'^'^{S^;M.'^), k > 0, be the constant speed parameterized 
unit circle, that is, 
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and let r < 2. Then the derivative of Ir['y] in direction a G C^^'^iS^;^^ with 

(T{ip) = a{ip)'y{ip)-^ + b{ip)j{Lp) 
is given by 

I'rb]i'^)i^) = / a{->p) dip - 2sm{-d)a{ip) = (x[-tf,^] * a)((^) - 2 sin(i?)a((/?) 



with 



'd := arccos ( 1 

2 



The proof of this leinina is postponed to the appendix. 

4. Proof of the Main Theorem 

Proof. Without loss of generahty we may assume that R = 1. Let V be the neigh- 
borhood of the unit circle defined in Theorem 3.3. The formula for /^ implies 
that for every 7 e V C (7*^+^(5^; M^) the mapping I^[y] is bounded as a map- 
ping from C''{S^;M.'^) to C'^(S'^;M). Thus, I'^[y] has a unique bounded extension 
Jrll]- C'^(5'^;K^) — )■ C'^(S'^;]R). Moreover, the mapping Jr is continuous with re- 
spect to the C'^+2-topology seen as a mapping from V to L{C''iS^;R^),C''{S^;R)) . 
In addition, we denote for 7 G V by Jr['y] the restriction of Jrl"/] to T^Ck- 

Denote now by 70 : 5^ — >■ M^ the constant speed parameterized unit circle. Define 
for given a G T^oCk the function Aa : 5^ — > M by 

Aa((^):=(a((^),7o((^)-') . 

Because 70 is a C°°-curve, it follows that Aa is C^ . Using Lemma 2.4 it follows 
that Acr(O) = and d^{A(T){Q) = 0. Denoting by C^{S'^;M.) the space of all C*^- 
functions a on the circle satisfying a(0) = a(0) = 0, it follows that A is a bounded 
linear mapping from T^^Ck to Cq{S^;M.). 

In addition, it follows from Lemma 2.4 that A is boundedly invertible with A^^ 
given by 

A^^a = a-JQ + 670 
with 

ip -1 p27T 



b{ip) = 5(0)(p + / a{T) dr and 5(0) = — / a{T) dr . 
Jo 27r Jq 

The expression for 6(0) is due to the periodicity of 5, which implies that 
= &(0) = 6(27r) = 27r6(0) + / a{T) dr . 



JQ 

Therefore A is in fact an isomorphism between T^gCk and Cq{S^;M.). 

According to Lemma 3.4, the mapping Jr[7o] evaluated at cr = aj^ + bjg G T^gCk 
can be written as 

Jr ho] (o-) =X[~^,i)]*a-2 sin(i?)a . 
Thus Jr [70] can be decomposed into 

Jr[7o] — B o to A , 
where the operator B: C''{S^;R) -^ C''{S^;R) is given by 

Ba = X[-'d,i)] * a — 2 sin(i?)a 
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and I is the embedding from Co(S'^;M) into (7*^(5'^; K). Lemma 6.1 (see Appendix) 
implies that the mapping a i— >■ X[-i5,i?] * a is compact and thus i? is a compact per- 
turbation of the identity. Therefore the Riesz-Schauder theory (see [17, Chap. X.5]) 
imphes that B has a closed range. 

Next we compute the kernel of -B. To that end we consider the mapping in the 
Fourier basis. A short calculation shows that in this basis the operator B is the 
diagonal operator that maps a sequence of (complex) Fourier coefficients (cj)jgz to 
the sequence (djCj)j^z, where 

'2(l-sin(i?)) ifi = 0, 



25MM_2sin(^) 



else. 



Because sin('i9) j^ 1 and sm{ji3) ^ j sm{d) whenever j e Z \ { — 1, 0, 1} (see Lemma 
6.2 in the Appendix), it follows that the kernel of B consists of the functions a of 
the form a{ip) — c_i exp(—iip) + ci exp(i(/?) for some c_i, Ci G C 

In the next step we show that the kernel of Jr [70] = B oio A is trivial. Therefore 
assume that a — c_i exp(— i-) + ci exp(?;-) e Cq{S^;M.) n KerS. Because a(0) = 0, 
it follows that c_i + ci =0; because a(0) = 0, it follows that — c^i + ci =0. 
Together, this shows that c_i = ci = 0, implying that the intersection of Ker_B 
with Cq{S^;M.) is trivial. Since A is an isomorphism this proves the injectivity of 

Mio]- 

We have thus shown that Jr[7o] = B o i o A is strongly closed and injective. 
Now note that the set of strongly closed and injective, bounded linear functionals 
between two Banach spaces X and Y is open with respect to the norm topology on 
L{X, Y). Because of the continuity of Jr and therefore Jr, this proves the existence 
of a neighborhood jq E U C Ck+2 such that Jr[j] is injective for every -f £ U. In 
particular, this proves the injectivity of Irlj] for every j € U. 

The local uniqueness follows directly from the Frechet differentiability of /^ and 
the injectivity of /^ on 14. D 

5. Conclusion 

In this article, we have shown a local uniqueness result for the circular inte- 
gral invariant on a neighborhood U of the circle. Note, however, that the derived 
result does not prove the local injectivity of the invariant on U. We still believe 
that also the injectivity holds, but the following difficulties arise using our line of 
argumentation: 

For an application of the inverse function theorem, we would require that the 
mapping Ir is Frechet differcntiable and its derivative /^ is an isomorphism of the 
corresponding tangent spaces. Although our results prove that /^ can be extended 
to an isomorphism Jr of C'^, we cannot use the inverse function theorem, as the 
mapping Ir is not Frechet differcntiable (and not even Gateaux differcntiable) from 
C'^ to C^ — in fact, our results do not even prove that /,, maps C*^ curves into C^ 
integral invariants. Conversely, seen as a mapping from C'^"''^ to C'^, the tangent 
mapping, though injective, is not a surjection near the circle. The main prob- 
lem regarding the differentiability of Ir is that our formula involves a composition 
mapping (see Lemma 3.1 and Theorem 3.3). 

Acknowledgment. The authors want to thank Peter Michor and Giinther Her- 
mann for their comments and suggestions, which helped to improve the article. We 
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6. Appendix 

Lemma 6.1. For every fc > 0, the mapping a i— >■ K§a := X[-^,-d] * a is compact as 
a mapping from C''{S^;R) to C''{S^;R). 

Proof. Denoting by Bi the unit ball in C'''('S'^;M), we have to show that the image 
of Bi under K^ is precompact in C'^{S^;M.). Applying the Arzela-Ascoli Theorem 
(see [17, Chap. III. 3]), we have to show that K^{Bi) is bounded and the first k 
derivatives of the functions in K§{Bi) are equicontinuous. The boundedness of 
K^{Bi) is obvious, K^ being a bounded linear mapping (of norm 2^). Now assume 
that a e Bi. Then ^j (Vi,a)(v5) = a^^''^\Lp+'d)-a'~^'^\ip-d). Because ||a(j)||oo < 1 
for all 1 < J < fc, it follows that d^p{K^a) is Lipschitz continuous with Lipschitz 
constant at most 2. Hence K^{Bi) is a precompact set. D 

Lemma 6.2. Lef < i9 < tt and j e Z \ {-1,0, 1}. Then s\n{jd) ^ j a\n{d) . 

Proof. Assume first that < i? < 7r/2. We show that in this case the equation 
sin(s) ~ ssin(i?)/?9 has the only solutions s = and s — ±-3. First note that the 
strict concavity of the sine function on the interval [0, tt] implies that on this interval 
we only have two solutions, namely and d. Moreover, the concavity of the sine 
implies that sin(?9)/i? > sin(7r/2)/(7r/2) = 2/tt, and therefore 7rsin(t9)/'!9 > 2. This, 
however, implies that the equation sin(s) = s sin(i?)/'i9 cannot have any solutions for 
s > TT, as the right hand side is strictly larger than 2. The fact that —d is the only 
negative solution follows by symmetry. In particular, setting s = jd, this proves 
the assertion in the case < d < -k 12. 

Now assume that 7r/2 < ^ < tt and let tp :— tt — ■&. Then 

sin(i?) = sin(7r — ip) — sin(V') . 

Now, if j is odd, then 

sin(j^) = sin(j7r — jip) = sin(7r — ^"0) = shi{jijj) . 

Thus j sin(i?) = sin(j'!9), if and only if jsin(-0) — sm{jip). Because < ip < 7r/2, 
the first part of the proof can be applied, showing that jsin(i?) ^ sm(ji3) unless 
j=±l. 

On the other hand, if j is even, we have 

sin(j'^) = sin(J7r — jip) ~ sin(— j?/') — — sm{jijj) . 

Thus, J sin(i?) — sin(j'i9), if and only if j'sin(?/;) = —shi(jijj). Now note that 
the equation — sin(s) = ssin^ip) /ip has only the trivial solution s = 0, because 
sm{ijj)/ip > 2/tt and the left hand side is negative for < s < tt. As a conse- 
quence, the equation j sin(V^) = — sin(j'(/;) only holds for j = 0, which concludes the 
proof. D 

6.1. Proof of Lemma 3.4. 

Proof. Let 7 e C'^+^(S'^; M^) be the constant speed parameterized unit circle. Then 
7(93) = {cos{ip),sm{ip)) , jiip) = (- sin(v3), cos{ip)) , 

7((^)-L = (sin((^), - cos((/3)) , jiip)-^ = (cos(^), sm{ip)) . 
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Because 7 is the unit circle, there exists ^ S S"^ such that 

P^{lp) ^ ip + -d , my{(p)^ip-d. 

Obviously, all assumptions of Lemma 3.3 are satisfied. It remains to calculate all 
the terms that appear in the expression of the Frechet derivative in Lemma 3.3 for 
the special case of the unit circle. In particular, we obtain 

^' = 117(^7) -7ll' = 2(1 -cosW), 



{i{Pi)n{Pi 

{i,i{Pi 
ii^nip-r 

{i{Pi)^,l{'m^ 

(7^,7(m^ 

(7(m^),7(p^ 

(7(rnT,)^,7(m^ 

(7(rnT,),7(m^ 

(7,7(to^ 

(7(m^)-L,7(p^ 

(7, 7(^7) -7 



-7) =-2cos(^)(l-cos(i?)), 

-7) =sin(i9)(l-2cos(t?)) , 

- 7) = 1 — cos(??) , 

- 7) == sin(i9) , 

- 7) = sin(?9) , 

- 7) = cos(i?) — 1 , 

- 7) = cos2(i?) - cos(i?) - sin^(i?) , 

- 7) = cos(^) — 1 , 

-7) = -sin(i?)(l-2cos(i?)), 

- 7) = 1 — cos(i9) , 

- 7> = - sin(?9) , 

- 7) = - sin(i?) , 

- 7) = cos2(i9) - cos(7?) - S\T?{d) , 
m^)) =2sin(i9) . 



In the following we treat the derivative in direction 07^ and bj separately. For the 
derivative in normal direction 07"*" we get 



P-y 



2/;[7](a7^)=2 / a^d^- 0(7^,7^7)^ -7(^7)^) 



4Pt){l{Pt) - 1, Up-,)) + a(TO^)(7(m^) ~ 7, 7(m^)) 



+ (7(^7) -7, 7(^7) > 



_L, (07-^ - a{p^)'j{pj)-^,j{p^) - 7) 



(7(P7)'7(P7)-7) 



(7(rn^)-7,7(m^) ) 



j_ (a-y-L - a(m^)7(m^)-L, 7(m^) ~ 7) 



(7(^7)' 7(w^)- 7) 



-(7(P7)'7("i7)-7) 



\/r^ - iliP-r) ~ 7, 7(^7) - 7)^ 

1^ (07-^ -a(p-y)7(p^)-L,7(p-y)-7) ^.^ 
A (7(P7),7b7)-7) ^^^^-^"^'--^ 

+ (a(P7)7(P7)^ - a-i^ilim^) - 7> + (7(^7) " 7, 0(771^)7(771^)-^ - 07^) 
(07-^ -a(m^)7(7n,^)-L,7(m^) -7) 



(7(TO^),7(7n.^)-7) 



-(7(^7)- 7,7(^17)) 
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Inserting the expressions for 7 equal to the unit circle that have been calculated 
previously, we obtain 

2/;[7](a7^) 

= 2 / a{ip)di}} 

(cos(i?) — 1)^ 
- sin(^) (a(p^) + 2a + a{m^)) ^TT\ ("(^7) + 2a + a{my)) 

(a + a(p^)) (cos(?9) - 1) (1 - 2 cos(i?)) - 2a(cos(i?) - 1) 



sin(^) 

+ {a{py) + a{m^)) (cos2(?9) - cos(i?) - sin^(^)) 

+ (a + a{m^)) (cos(i9) - l) (l - 2 cos(i9))) 



= 2 / a(V')dV 

cos(??) - 1 



(2a(p^) + 4a + 2a(m^) - 2a(p^) + 4cos(^)a - 2a(m^)) 



sin(^) 
Jip-tf sin(i?) 



= 2 / a(-0)dV-4asin(??) . 

Inserting the formulas above in the expression for /^ [7] (67) yields 
2/;[7](67) = 2 P 6(^)(7(V'),7W^)rfV' 

- {bJniPi)^ - liiT^i)^) + iliPi) - 7,^(^7)7(^7)^) 

- (7(^7) -7,^-7)7(^7)-) + (7(P7) - 7, 7(P7)") ^'' '( • (;;)l!f(;!^7!';) ' '^ 



U\ (^')' - K'7i7)7(m^), 7(n^7) " 7) 
'-'-'^''"'-^^ > (7(m,),7(m,)-7) 



(7(m^)-7,7(m^)-^)- 



^r4 - (7(p^) - 7, 7(m^) - 7)2 

-(7(P7),7('™7) -7) 
+ {KPi)i{Pi) - &7,7(nT-7) - 7) + (7(^7) - 7,^(n^7)7(»7^7) - 67) 

-(7(P7)-7,7('™7))j 
= - + (1 - cos(i?)) (6(p^) - b{m^) + h- b{p^) -b + b{m^)) 

- f(l - 2 cos(z9)) (5 - b{py) + b{py) - b{m-y) + b{m^) - b) + {b - bfj 
= 0. 

Therefore, 

/^[7](a7-^ + bj) = X[-i5,,5] * a -2sin{i3)a . D 



(^7-^(^7)7(^7), 7(^7) - 7) /,„ 
(7(lJ7),7(P7)-7) 



{bj - b{m^)j{Tn^), -f{m^) - 7) 
(7(m^),7(m^)-7) 
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